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^ I of its phase diagram. We find a dense spectrum consisting of an almost unaltered 

repetition of the glueball spectrum of the pure gauge theory, and additional bound 
J2 ' states of adjoint scalars. For the parameters chosen, the model represents the effective 



finite temperature theory for pure SU(2) gauge theory in four dimensions, obtained 



D , after perturbative dimensional reduction. Comparing with the spectrum of screening 

^ ' masses obtained in recent simulations of four- dimensional pure gauge theory at finite 

temperature, for the low lying states we find quantitative agreement between the full 
^ ' and the effective theory for temperatures as low as T = 2Tc. This establishes the 

■ model under study as the correct effective theory, and dimensional reduction as a 

viable tool for the description of thermodynamic properties. We furthermore compare 
the perturbative contribution ~ gT with the non-perturbative contributions ~ g^T 
and ~ g^T to the Debye mass. The latter turns out to be dominated by the scale g'^T, 
whereas higher order contributions are small corrections. 
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1 Introduction 



The thermodynamic properties of non-Abehan field theories at finite temperature are 
known to be non-perturbative due to the confining nature of the Matsubara zero mode 
sector. Lattice simulations of the full problem are numerically expensive and espe- 
cially difficult if fermions are present. A way to circumvent both problems is to per- 
turbatively integrate out the non-zero Matsubara and momentum modes of order ~ T 
and higher to arrive at a purely bosonic, three-dimensional effective theory describing 
the static quantities of the original theory [|l|. Three-dimensional gauge theories are 
non-perturbative, but readily accessible by lattice simulations. This programme has 
been successfully carried out in studies of the electroweak phase transition, where the 
phase structure and nature of the transition were computed from the dimensionally 
reduced theory, both by numerical [|] and analytical methods. The results have 
been found to agree to remarkable accuracy with the answers obtained recently from 
four-dimensional simulations [Q. 

In view of the experimental search for the quark gluon plasma in heavy ion collisions 
at RHIC and LHC it is of particular interest to see to what extent this programme 
can be carried over to QCD. The reduction step for a SU(N) gauge theory with Nf 
flavours of fermions has been carried out to two loops in (see also |^). The effective 
theory obtained in this case is a three-dimensional SU(N) adjoint Higgs theory, with 
the scalars corresponding to the electric gauge potential, Aq, in the unreduced theory. 
In the framework of dimensional reduction Nf only modifies the parameters of the 
effective, bosonic theory. To be specific, and to compare with existing four-dimensional 
work, we discuss in this study the simplest case of iVj- = 0, = 2. 

Simulations of the three-dimensional SU(2) adjoint Higgs model revealed that this 
theory has a phase diagram with continuously connected Higgs and confinement regions 
0, § , which are partially separated by a line of first order phase transitions. This is in 
contrast to the four-dimensional SU(2) Yang-Mills theory at finite T, which displays a 
second order phase transition ^j. Hence the deconfinement phase transition cannot be 
described by the effective theory. Based on the apparent convergence properties of the 
reduction step, however, it has been argued that the effective theory still yields the 
correct description of static correlation functions at temperatures sufficiently above the 
critical temperature, Tc, and hence valuable information about the plasma state. 

In recent publications ||^, |TD[ it was established that dimensional reduction indeed 



works in the sense that the spectrum measured in four- dimensional simulations at finite 
temperature can be classified by the symmetry group of a purely three-dimensional 
theory. This does not yet determine the precise form of the effective theory. Earlier 
simulations of the SU(2) adjoint Higgs model have concluded it to be a good effective 
theory by comparing the static potential |TT| and gauge-fixed propagators [|T^ between 
the full and the effective theory. The question of the dynamical role of the Aq and 
whether there is a further hierarchy gT <^ g^T in the reduced theory, however, is still 
not precisely answered. The purpose of the present paper is to quantitatively test 
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the agreement of a wide range of gauge-invariant correlation functions between the full 
and the effective theory at temperatures not much larger than the critical temperature. 
We compute correlation functions of the three-dimensional SU(2) adjoint Higgs model 
on the lattice and extract the corresponding screening masses for various quantum 
number channels. These are then compared with the four-dimensional calculation of 
finite temperature SU(2) Yang-Mills theory in |TD|. We find indeed that dimensional 



reduction works remarkably well, even at the quantitative level. Furthermore, we are 
able to identify the lightest gauge-invariant state to consist of Ao, which hence may 
not be integrated out. 

The continuum action of the SU(2) adjoint Higgs model is given by 

S = J d^xi^^Ti {Fi,F,j) + Tr (Di^D,^) + m^Tr {^^) + X-^iTi (w))'} , (1) 

where Fij = diAj — djAi + ig[Ai, Aj], Diip = diip + ig^lAi, ip] and Fj^, Aj, and ip are all 
traceless 2x2 Hermitian matrices {(p = ip"'aa/2, etc.). The physical properties of the 
theory are fixed by the two dimensionless ratios 

93 93 

In the framework of dimensional reduction these parameters are completely determined 
by the four-dimensional gauge coupling, g"^, and the temperature, T. Since g"^ is a 
running coupling its value is in turn fixed by a renormalisation scale Aj^. The details 
of the derivation of x,y to two loops can be found in Here we merely quote 
the results to furnish the connection between the three-dimensional parameters and 
the four-dimensional, finite T situation. Choosing the renormalisation scale as in |^ 
and expressing through the critical temperature as measured on the lattice 0, 
Tc = 1.23(11)A^, one has 

2 _ 10.7668 _ 0.3636 r ^ - ^ ^ 

^3"ln(8.3T/T,)^' ln(6.6T/Tj' ^^"^^ " 9^ + 4^^ + 

With these equations, specifying T/T^ completely fixes the parameters x,y of the re- 
duced model. 

In this paper we study specifically the temperatures T = 2Tc, T = 4Tc and T = 5Tc 
(with the corresponding values of x and y given in Table [^). Note, however, that these 
parameters lie in the Higgs phase of the model, whereas the connection to four dimen- 
sional physics is only valid in the confinement phase |T^, The dimensional reduction 
programme would be futile were it not that the confinement phase is metastable in this 
region of the phase diagram. Starting with a configuration in the confinement phase, 
we may re-thermalise and make many measurements at the appropriate parameters 
before a tunnelling transition occurs to the Higgs phase. In practice, the phases are so 
strongly separated that no tunnelling occurred in our simulations. 
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In order to discuss correlation functions of the effective theory, it is most convenient 
to ignore the finite T context and to discuss the properties of the 2 + 1 dimensional 
Higgs model by itself, as we shall do in the following two sections. This will also enable 
us to compare the dynamics of the SU(2) adjoint Higgs model in its confining phase 
with that of the 2 + 1 dimensional pure gauge theory and fundamental Higgs model, 
whose spectra are known in some detail IT^ O] , and thus to learn about some general 



features of three-dimensional confinement. In Sec. 2 we introduce the corresponding 
lattice action, the operators considered and some technical details about our numerical 
simulation. In Sec. 3 we present our detailed numerical results for the spectrum of the 
model, including checks for finite size effects and a numerical continuum extrapolation. 
Then we return in Sec. 4 to an interpretation and discussion of our results in the context 
of finite temperature pure gauge theory, before giving our conclusions. 

2 Lattice methods 

2.1 Action and parameters 

The discretised form of the action ([l|), with leading order lattice spacing corrections of 
0{a), is 

S = /?E (l-^Trf/,,(x))+2 5:Tr((p(xMx)) 

x,i>j \ ^ / X 

- 2KY,Ti{ip{x)U^ixMx + fia)U\x)) + Xj2i'2Ti{if{xMx))-lf. (4) 

(The lattice scalar field has been rescaled relative to the continuum.) The parameters 
of the continuum and lattice theory are up to two loops related by a set of equations 
specifying the lines of constant physics in the lattice parameter space {/?, k. A} |[T6[] , 
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(5) 



where S = 3.17591. For a given pair of continuum parameters x,y these equations 
determine the lattice parameters k, A as a function of lattice spacing, and hence govern 
the approach to the continuum limit, /? — >• cx). It is along such trajectories that 
we should observe a scaling of dimensionful quantities such as the masses. Due to 
superrenormalisability of the theory, the above perturbative relation is exact in the 
continuum limit. Based on previous experience M we expect it to be accurate for all 
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P > 6. As a consequence of the presence of scalar fields in the action, the leading order 
lattice spacing corrections to mass ratios should be linear in a. 



2.2 Blocking of fields 

Since we work in the confinement region of the theory, the physical states are expected 
to be bound states, which are spatially extended. In order to improve the projection 
of any operator onto such states, it is necessary to smear the operators in the spatial 
plane so that their extension becomes comparable to the size of the states. This can 
be achieved by using "blocked" field variables instead of the elementary ones in the 
expressions for the operators. We use blocked link variables at blocking level n as 



proposed in [|T7 



For the scalars, we adapt the procedure employed in to the case of adjoint fields. 



Note that we will take correlations in the 3 direction and so our blocking always remains 
within the (1, 2)-plane to avoid spoiling the transfer matrix of the theory and the 
positivity properties of the correlation functions. In addition, we consider blocked 
scalar fields Lp^''^''-'\x) with non-local contributions from one spatial direction j only: 

<^(-.i)(x) = ^((^("-^'^■)(a;) + f/f-'^(x)(^("-i'^'nx+j)f/J""'^^(a;)) (8) 

+ uf'~'^\x-])^'^-''^\x-])uf'-'\x-])], J = 1,2. 



2.3 Basic operators 

The operators we use in our simulations are constructed from several basic gauge- 
invariant operator types, which we now list before specifying the quantum numbers. 
We have operators involving only scalar fields or products of scalar and gauge field 
variables, 

R{x) = Tt{^\x)) 
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Liix) 



A-.(x) 



-- Tr {Ui{x)Ui{x + + 2i)Ul{x + i)Uj {x)Uj{xMx + ])Uj{x] 
-- Tr {Ui{x)Ui{x + i)<^{x + 2i)f//(x + i)U^ {x)Uj{x)i^{x + ])Uj{x)<^{x] 
Other mixed operators are 



(9) 



= Ti{ip{x)Uij{x)), 
B^{x) = TtT" lUijix) - ip{x 



Tr {(p{x)Uij{x)) 



(10) 



The naming of the operators becomes apparent when considering the continuum hmit, 
in which B\\{x) ~ Tt {(p{x)Fij{x)) gives the projection of the field strength along the 
scalar field and B±{x) is the projection perpendicular to it. 

Further, we have loop operators constructed from link variables only, 



a 



1x1 



(x) = Tr {u,{x)U,{x + i)Uhx+])Uj{x)) , z,j = 1,2, z ^ j, (11) 



and in addition to the elementary plaquette C^^^, we also consider squares of size 2x2 
as well as rectangles of size 1x2, 1x3, 2x3. Hence, we have five versions of operators 
consisting of closed loops of gauge fields, viz. 



f^lxl ^2x2 ^1x2 ^1x3 ^2x3 



(12) 



Another pure gauge operator useful to probe the confining properties of the theory is 
the Polyakov loop along a spatial direction j, 



p^^Hx) 



L-l 

Tr n Uj{x 

m=0 



J 



1,2. 



(13) 



Apart from its relevance for confinement, the Polyakov loop also plays an important 
role in the study of finite size effects. Although single torelons arising from correlations 
of the Polyakov loop Pj^^ can neither contribute directly to correlations of operators 
C^^^, . . . ,C'^^^, nor to our mixed operators, they may well couple indirectly on the 
grounds that they share the quantum numbers with the infinite volume states. Fur- 
thermore, torelon-antitorelon pairs are known to give rise to finite-volume effects in 
glueball calculations |jl8|. Apart from studying correlations of Polyakov loops per se, 
we have thus constructed single-torelon and torelon-pair operators in the 0^ and 2^ 
channels and studied their correlations as a safeguard against finite-size effects. 



2.4 Quantum numbers 

Operators with definite quantum number assignments may be constructed from the 
above operator types by taking linear combinations with appropriate transformation 
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properties, i.e. by applying projection operators for the various irreducible represent- 
ations of the symmetry group of the lattice slice upon which the operators are con- 
structed. 

In a three-dimensional Euclidian gauge theory, the slice is a plane with an S0{2) 
rotational symmetry and discrete Z2 parity reflection and charge conjugation symme- 
tries. Rotations are generated by R{6) under which operators of angular momentum 
j gain a phase exp[ij^]. Parity in two spatial dimensions reflects operators in one 
coordinate axis; P : {x,y) —>■ {x, —y); reflections in other axes are related to this by 
rotations. Charge conjugation acts as complex conjugation on the matrix valued fields, 
C : Ui{x) U*{x),ip{x) ^*{,x). Physical states and operators are characterised by 
their eigenvalues under these operations, the quantum numbers J^^ . In contrast to 
the fundamental Higgs model, there are no C = —1 states in the SU(2) pure gauge 
theory or the SU(2) adjoint Higgs model. In this model SU(2) charge conjugation is 
equivalent to a global gauge transformation, and hence any gauge-invariant operator 
is even under this operation. We therefore omit the index C from now on, implying 
that all states have C = +1. 

With an adjoint scalar field present, there is an additional Z2 symmetry correspond- 
ing to reflections Lp — ^ —ip, whose eigenvalue we call R. Clearly operators containing 
an even number of scalar fields, such as A^, will couple to i? = +1 states, and those 
with odd, such as A^^-, to i? = — 1. We thus classify our operators and the states they 
couple to by J^. 

On the lattice the symmetry group is broken from SO{2)®Z2{P)®Z2{R) to the point 
group ® Z2{R), which restricts the allowed rotations to R{On) where 9n = mi/2. 
Strictly speaking, we should thus classify our states by the irreducible representations 
of rather than by J^. Since we are really interested in continuum physics and, as we 
shall see, our data reproduce the continuum symmetries within the statistical errors, 
we prefer to keep the continuum notation. 

The fact that the two-dimensional parity operator reflects only one spatial direction, 
and hence the angular momentum operator changes sign under P, has an important 
consequence for the spectrum: in 2 -|- 1 dimensions all states with J 7^ come in 
degenerate pairs of opposite parity (see, e.g., [0). This statement of "parity doubling" 
is based on the continuum rotation group, whereas a square lattice with a finite volume 
only admits rotations in units of 7r/2. Thus, although we should continue to observe 
parity doubling of 1^ and 1^ on the lattice, the doubling in the J = 2 sector is 
lost, and will only be recovered in the infinite volume limit. Comparing the masses 
of the 2^ and 2^ states is thus a useful measure of whether our lattice is free from 
discretisation and finite volume effects. In this study all simulations were done on 
lattices large enough that parity doubling was, within statistical errors, in accordance 
with continuum, infinite volume expectations. 

An interesting feature of the SU(2) fundamental Higgs model is the near replication 
of the pure gauge theory glueballs within its spectrum. Whereas one might expect the 
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excitations of the former to be a mixture of gluonic and scalar degrees of freedom, it 
is found instead that within the spectrum there are some essentially purely gluonic 
states. Given their remarkable decoupling from the scalar fields, it is then unsurprising 
that the masses of these match that of the pure gauge model almost exactly |TB[ . One 



of our main interests is to study whether the glueball spectrum survives the addition 
of adjoint scalar fields similarly unchanged. From the known spectra of the pure gauge 



theory [0] and of the fundamental representation Higgs model we expect the spin 
and spin 2 channels to contain the lowest glueball masses, so in order to carry out the 
necessary mixing analysis, we have constructed large bases of different operator types 
for these two channels: 

0^ channel: 
R: R{x) 
L: Li{x) + L2{x) 



B: B^{x) 

C: symmetric combinations of C^^\ C^^^, C^^^, C^^^, C^^^ 

P,: Pi'^^x) ■ Pi'^^x) 

T: {Pl'\x)y + {Pi'\x)f 

2^ channel: 
R 



L 

A: 
C 
P 
T 



pL"^ (x) = 1{Tt b^"'^) (x)] - Tr (x)] } , n > 1 

Li(x) - L2{x) 

ELi e^^'^Pl^n) {A^x) + PAUx) 



antisymmetric combinations of C^^^, C^^^, C"^^^ 
P[^\x) - P!^^\x) 
(P['\x))' - (P^'\x)\ 



In these expressions R{6n) rotates the following brackets by 6n around x, and P denotes 
the parity operation. The labels P, Pd refer to single torelon operators, whereas T 
couples to torelon pairs in the respective channels. 

In principle, similarly extended bases are possible for the other channels as well. 
This requires a lot of memory and computer time, and furthermore the glueballs in the 
other channels are anticipated to be rather heavy ||15|, Q , and a mixing analysis would 
be difficult. The dynamics of mixing is moreover expected to be independent of the 
precise quantum number channel so it should be sufficient to study the 0^ and 2+ 
as detailed examples. For the other channels we therefore limit ourselves to just one or 
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two operator types B\\ and/or A. We only give two more explicit examples; the other 
combinations follow straightforwardly. 

Ol channel: 

A: Y.n=iR{On)[A'{^{x)-PAUx)) 
ll channel : 



For all channels we construct operator versions from elementary as well as blocked 
fields at various blocking levels. We then take zero momentum sums by averaging over 
a timeshce, (f){t) = l/L'^J2xi,x2't>{x)^ where (t>{x) denotes a generic operator at some 
blocking level. 

2.5 Matrix correlators 

In order to compute the excitation spectra of states with given quantum numbers, 
we construct matrix correlators by measuring all cross correlations between different 
types of operators at several blocking levels. The correlation matrix can be diagonalised 
numerically following a variational method. For a given set of operators {(pi} we find 
the linear combination that minimises the energy, corresponding to the lightest state. 
The first excitation can be found by applying the same procedure to the subspace 
which is orthogonal to the ground state. This may be continued to higher states so 
that we end up with a set of N eigenstates i = 1, . . . , N} given by 

AT 

= E^ifc'^fcW- (14) 

fc=i 

The coefficients a^fc quantify the overlap of each individual operator 0^ used in the 
simulation onto a particular, approximate mass eigenstate For a complete basis of 
operators this procedure is exact. In practice, the quality of the approximation clearly 
depends on the number N of original operators and their projection properties. The 
eigenstates already determined are removed from the basis for the higher excitations, 
so that the basis for the latter gets smaller and the corresponding higher states are 
determined less reliably. A more detailed discussion of the calculation and of checks 



of its stability can be found in |15|, where the same procedure was applied to the 
fundamental Higgs model. 

In order to distinguish glueballs and scalar bound states, the most elaborate calcula- 
tion is done in the Ot and 2t channels. For Ot our basis consists of = 32 operators 
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0fe G {R,L,A,B,C,P,T,Pd}, for 2+ we have = 31 and 0^ G {i?, L, A, C, P, T}, 
each operator considered for at least two blocking levels. This typically enables us to 
extract the four or five lowest lying states in those channels reliably. In the Ol, 2l, 
ll, ll channels we work with = 9 operators 0^ G {A,B}, and in the remaining 
channels we have = 4 and 0^ ^ {^}- Correspondingly, in the smaller bases we can 
only extract the ground state and, in some cases, the first excitation reliably. Finally, 
in order to compute the string tension we have also diagonalised a separate basis of 
Polyakov loop operators at = 4 different blocking levels. 

2.6 Simulation and analysis 

The Monte Carlo simulation of the lattice action, Eq. (^, is performed using the same 
algorithm as in previous work 0. The update of gauge variables employs a combination 
of the standard heatbath and over-relaxation algorithms for SU(2) |]19|, whereas the 
scalars are updated by a suitable adaption of the algorithm for fundamental scalars 
suggested in ||2^ and employed in . In order to account for the quadratic dependence 
of the hopping term on the link variables, the link update for the pure gauge action 
has to be supplemented by a Metropolis step. 

We define a "compound" sweep to consist of a combination of one heatbath and sev- 
eral over-relaxation updates of the gauge and scalar fields. In this study we performed 
5 over-relaxation steps for every heatbath update. Measurements are taken after every 
such compound sweep. We have gathered between 5 000 and 30 000 measurements, 
depending on the lattice sizes which range from 24'^ for the coarsest lattice to 48^ on 
the finest one. 

All our mass estimates are obtained from measured correlation functions of operators 
$j in the diagonalised bases defined in the preceding subsection. Correlated fits over 
some interval [^1,^2] are performed to the ansatz for the asymptotic behaviour of the 
correlation function on a finite lattice, for large T, 

(<l'l(t)$i(0)) = Ai (e-"^'»* + e-"^'^^-*)) , (15) 

where i labels the operator, and T denotes the extent of the lattice in the time direction. 
In cases where the effective masses did not show plateaux long enough to perform a 
correlated fit, we used the effective masses for our mass estimates. 

Statistical errors are estimated using a jackknife procedure for which the individual 
measurements have been accumulated in bins of 100 — 250 sweeps. 

3 Numerical results 

We have performed simulations at temperatures T = 2Tc, T = 4Tc and T = 5Tc 
corresponding to three points {x, y} in the metastable confinement phase of the model, 
as detailed in Table 
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T 




y 


/3 = 


9 




/? = 


12 




(3 = 18 




0.141 


0.185 




T = 28^, 


38^ 


L^- 


T = 


38^ 






0.111 


0.228 








L2. 


T = 


38^ 




5T, 


0.104 


0.242 


L2. 


T = 243, 


343 


L2. 


T = 


343 


L^-T = 483 



Table 1: The lattice parameters and sizes simulated. 



T = 2Te 


(3 = 9 


(3 = 12 




L = 38 


L = 28 


L = 38 


aMp{L) 


0.849 (13) 
0.151 (1) 
0.339 (3) 


0.627 (6) 
0.152 (1) 
0.341 (2) 


0.462 (4) 
0.112 (1) 
0.336 (2) 



Table 2: Polyakov loop masses and string tensions atT = 2Tc. 



The most elaborate test of finite volume effects as well as an explicit continuum 
extrapolation for the lowest states are done for T = 5Tc. Experience of the pure gauge 
model |T^ suggests that to avoid contamination of the spectrum by light torelonic 
states the lattice size should be L > 24 at /3 = 9. We thus simulate three different 
lattice spacings (the third one being half of the first) while approximately maintaining 
this physical volume {i.e. L/ (3 = 24/9, using the lattice spacing as defined by the bare 
coupling). We find all these couplings to lie in the scaling window of the theory which 
permits a continuum extrapolation of the low lying masses. To verify that finite volume 
effects are indeed small we also simulate a second, larger volume at the coarsest lattice 
spacing. 

At T = 2T(. we simulate two lattice spacings and physical volumes, and at T = 4Tc 
just one as we are by now confident that lattice effects are small. In these cases 
continuum extrapolation is not possible and we thus use the masses measured at the 



T = 4T, 


(3 = 12 




L = 38 


aMp{L) 


0.476 (4) 
0.114 (1) 
0.342 (2) 



Table 3: Polyakov loop masses and string tensions atT = AT^. 
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T = 5T, 


(3 = 9 


(3=12 


(3=18 




L = 34 


L = 24 


L = 34 


L = 48 


aMp{L) 


0.80 (1) 
0.155 (1) 
0.349 (2) 


0.554 (5) 
0.155 (1) 
0.349 (2) 


0.426 (4) 
0.114 (1) 
0.342 (3) 


0.257 (6) 
0.075 (1) 
0.338 (4) 



Table 4: Polyakov loop masses and string tension atT = hT^ 



smallest lattice spacing as approximations of the continuum result. 

In order to perform continuum extrapolations it is necessary to have some definition 
of the lattice spacing with which to scale the dimensionless lattice masses. The simplest 
measure comes from the bare gauge coupling, a = A/{(3gl). Such a formula is liable 
to quantum corrections, which in some cases translate as large 0{a) corrections to the 
naive scaling of the masses A non-perturbative definition of the lattice spacing 



should in principle avoid this, and in pure gauge theory it is customary to scale all 
lattice masses in units of the string tension. We thus discuss this quantity first. 



3.1 The string tension 

As explained previously, a Polyakov loop in a spatial plane couples to a fiux loop 
state (or torelon) that winds around the periodic boundaries of the finite volume. The 
exponential fall-off of its correlator is related to the mass of the flux loop, according 
to 

^ (Pf )(x)Pf )^(0)) ~ e-«^^(^)*, aMp{L) = a^aLL. (16) 

X 

Such a flux loop state can be easily identifled through its energy scaling linearly with the 
size of the lattice, as seen in the entries labelled by P in Tables ^ ^ and In contrast to 
the fundamental Higgs model, there are no fundamental representation matter fields in 
the action to screen the colour fiux of static charges leading to a breaking of fiux tubes 
longer than a screening length |2^ . On the contrary, a fiux tube between fundamental 
charges will persist to infinite separation as in the pure gauge theory, and a string 
tension can be defined in precisely the same way by the slope of the static potential at 
infinite separation^. Accordingly, periodic fiux loops winding through the boundaries 
of the lattice will not be screened, and the coefficient (Jl corresponds to the string 
tension on a finite volume. An estimate for the string tension in infinite volume is then 
provided by the relation 

aVoo = ciVl + (17) 

^ To avoid confusion, we note that this is the string tension in the 2 + 1 dimensional theory, 
corresponding to the spatial string tension in the 3 + 1 dimensional theory at finite temperature. 
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Figure 1: The coefficients aik of the operators used in the simulation for the five lowest 
0^ eigenstates for 28 operators at (3 = 9 (left), and for 32 operators at j3 = 18 (right). 



We have diagonalised a basis of Polyakov loop operators and the results are summarised 
in Tables As we shall see with the glueball masses, the string tension is barely 
affected by the presence of the scalar fields. It also shows little dependence on the 
scalar parameters x, y. 

Expressing the string tension in units of the bare coupling, we may (linearly) ex- 
trapolate the results at T = 5Tc to the continuum, where we find ^/Coo/gi = 0.326(7), 
which compares to ^Jo^j g\ = 0.3353(18) in pure gauge theory [|14|. That the slope 
of this extrapolation is small indicates that the difference between using a tree level 
definition of the lattice spacing, and a non-perturbative one is slight, and in most 
cases any improvements in using the string tension as opposed to the bare coupling for 
setting the scale are outweighed by the additional statistical errors on the former. 

For this reason, in the following tables we scale and quote masses in units of the bare 
continuum coupling to allow rapid comparisons of measurements at different f3. 

3.2 The mass spectrum 

A detailed account of the mass estimates for the lowest states computed on the various 
lattices is given in Tables |5HlO|, together with the dominant operator content of a given 



mass eigenstate for the channels with the largest operator bases. 

Let us discuss in detail only the T = 5Tc parameter set, the procedure being en- 
tirely analogous for the others. As an example for the identification of the operator 
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T = 2Tc 


P = 9 


/? = 12 


state 


L^-T = 383 


L^-T = 283 


.T = 38^ 


ot 


M/5(| Ops. 


M/5(| Ops. 


M/^l Ops. 


$3 

$4 
$5 
$6 
$7 

$8 


0.747 (7) R,L,A 
1.629 (27) C, B± 
2.007 (45) A,{R,L) 
1.913 (45) P 
2.22 (Q) C,B± 
2.80 (13) A 
2.94(17) C, (5^) 
3.11 (48) (P) 


0.743 (7) R,L,A 
1.411 (18) P 
1.613(23) C,Px 
1.980 (36) A,{R,L) 
2.36(6) C,Px 
2.47 (9) P, 
2.75 (11) A 
2.86(13) C,B^ 


0.744 (9) P,L,A 
1.416 (15) P 
1.623 (21) C,P^ 
1.998 (27) A,{R,L) 
2.34 (4) C, (Px) 
2.38 (5) Pd 
2.79(7) C,B^ 
2.83 (7) (A) 


2+ 


M/(73" Ops. 


M/gl Ops. 


M/(73^ Ops. 


$1 

$2 
$3 

$4 
$5 
$6 


1.94 (5) P 
1.92 (5) P 
2.27(6) P,L,A 
2.61 (11) C 
2.60 (29) (P) 
3.15 (19) C 


1.440 (22) P 
1.454 (18) P 
2.24(5) P,L,A 
2.34 (8) (P) 
2.70 (12) C 
3.05(15) C,T 


1.401 (15) P 
1.419 (18) P 
2.30(4) P,L,A 
2.37 (5) (P) 
2.65 (5) C 
2.98(8) C,T 



Table 5: Mass estimates and dominant operator contributions in the 0^ and 2^ chan- 
nels at T = 2Tc. The dominant operator types contributing are denoted with (without) 
parentheses if a-ik < (>)0.5. 



content of a mass eigenstate, we display the coefficients aik {cf. Eq. (|D) of the five 
lowest mass eigenstates in the 0|-channel in Fig. |l] for our largest and smallest lattice 
spacing. The operator contributions listed in Table ^ are those which have coefficients 
Ojfc larger than 0.5. If the state has weak overlap with all operators of our basis and the 
dominant coefficients are smaller than 0.5, we indicate this by parentheses in the table. 
The ordering of states $j was obtained during the diagonalisation procedure from the 
effective mass of the corresponding correlator on the ffist timeslice. 

we 



As in the confinement phase of the fundamental representation Higgs model |15 



find here a dense spectrum of bound states also. In it was found that the Yang- 
Mills glueball spectrum is repeated almost identically in the model with fundamental 
scalar fields, with additional bound states of scalars, and little mixing between them. 
An interesting question is whether the situation is similar in the case with adjoint 
scalar fields. Fig. |l| (left) shows that this is indeed so. The ground state is clearly 
dominated by scalar operators, whereas the ffist as well as the fourth excited state 
have almost exclusively gluonic contributions and hence are identified as glueballs. 
Fig. |I| (right) is the same for our finest lattice, but now with four blocked versions 
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1 — 11c 


J 

/ ; - 


- y 


-1 1 

/ ' — J- - 


state 




Li — ZO 


Ij — OO 


Ot 






3 24 fl2) 


Ol 


1.229 (14) 


1.220 (14) 


1.118 (12) 




2.14 (6) 


2.14 (6) 


2.16 (4) 








3.04 (12) 




3.00 (5) 


3.01 (5) 


2.93 (11) 


11 


2.02 (2) 


2.01 (2) 


1.99 (3) 


ii 


1.95 (5) 


1.87 (8) 


1.92 (5) 


2+ 


2.41 (9) 


2.57 (9) 


2.52 (6) 


K 


2.35 (3) 


2.36 (6) 


2.31 (5) 



Table 6: Masses M/ for the other channels atT = 2Tc. 



T = 4Te 


P = 


12 


state 


L'-T 


= 38^ 


o| 


M/gl 


Ops. 


$1 


0.945 (9) 


R,L,A 


$2 


1.425 (15) 


P 


$3 


1.626 (18) 


C,B^ 


$4 


2.151 (30) 


A,{R,L) 


$5 


2.33 (4) 


C,B^ 


$6 


2.33 (4) 


(P) 


$7 


2.43 (5) 


Pd 


$8 


2.84 (7) 


C 


n 




Ops. 




1.428 (21) 


P 




1.452 (18) 


P 


$3 


2.42 (4) 


R,L,A 


$4 


2.34 (5) 


(P) 


<f5 


2.63 (6) 


C 


$6 


3.02 (8) 


C,T 



Table 7: Mass estimates and dominant operator contributions in the Oljl and 2^ chan- 
nels atT = 4Tc. The dominant operator types contributing are denoted with (without) 
parentheses if ttik < (>)0.5 
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T = 4Te 


(3=12 


state 


L = 38 


o; 


3.48 (12) 


Ol 


1.323 (12) 




2.30 (4) 


n 


3.47 (24) 


K 


3.05 (11) 


11 


2.00 (3) 


1- 


2.05 (4) 


2+ 


2.57 (5) 


2; 


2.45 (5) 



Table 8: Masses M/g^ for the other channels at T = AT^. 

of the operator i?^ (c/. Eq. (|T0|)) included in addition to the operators used for the 
left plot. At first sight this appears to violate the "non- mixing" between scalar and 
gluonic states as observed with the smaller basis. Upon expanding and taking the 
trace in Eq. ([T0|) , however, one readily sees that the operator is in fact a sum of terms 
of the form B± ~ Tr (f/^) + 5^ + Tr (<^^)i?^. Of these, the first term mixes with our 
other loop operators and has projection onto the glueball states since (Tr (f/^) has the 
same representation content as Tr(f/)), whereas the other terms project onto higher 
excitations (they are of higher dimension) and thus die out asymptotically. Hence, 
the operator asymptotically behaves as a C-type operator which mainly projects 
onto glueballs, but very little onto scalar states, as is corroborated in Fig. ||. This 
figure further shows that the decoupling of the gluonic sector is stable under variations 
of the lattice spacing by a factor of two. As we shall see, with the lattice spacings 
considered we are quite close to the continuum, and the observed effect should persist 
in the continuum limit. This picture is repeated in the 2^ channel. 

The spectrum of higher excitations becomes increasingly dense with increasing en- 
ergy. For the 0^ and 2+ channels, the eigenstates beyond $6 become increasingly 
difficult to identify, and their total overlap with the operators in the basis becomes 
smaller. Further, they are closer in energy and within our statistical errors appear 
often degenerate. An unambiguous resolution of these states would require still larger 
bases as well as more statistics. We have therefore chosen not to quote them as reliable 
results. 



Another word of caution concerns the states 0_,_, 1^ and 1_^. As described in Sec. ^ 



the bases for these states consist exclusively of four blocked versions of the operator 
type A. It is conceivable that these operator sets have very poor projection onto some 
of the sought ground states, leading to an overestimate of their masses. Indeed, such 
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T = 5Tc 




/5 = 


9 




state 




= 343 


L^-T-- 


= 243 


ot 

+ 


Mlql 


Ops. 


Mlql 


Ops. 




1.010 (9) 


R,L,A 


1.010 (12) 


R,L,A 


$9 


1.61 (5) 


C 


1.24 (3) 


P 




1.82 (7) 


p 


1.64 (3) 


C, B± 


~±r 4 


2.19 (14) 


L,A 


2.16 (11) 


L,A 




2 29 fl6) 


C 


2 05 fl3) 


Pi 


u 


2.73 (15) 


(P) 


2.36 (18) 




2+ 




Ops. 


M/gl 


Ops. 


$1 


1.795 (32) 


P 


1.215 (25) 


P 


$2 


1.798 (32) 


P 


1.253 (25) 


P 


$3 


2.35 (8) 


R,L,A 


2.25 (10) 


(P) 


$4 


2.68 (11) 


C 


2.45 (15) 


R,L,A 


*5 


2.67 (11) 


{P) 


2.57 (11) 


C,T 


$6 


2.87 (18) 


C 


2.86 (11) 


C,T 


$7 


3.34 (24) 


R,L,A 


3.17 (17) 


R,L,A 




T = 5Tc 




12 


P = 


18 


state 




= 343 


L^-T-- 


= 483 


oi 


M/gl 


Ops. 


M/gl 


Ops. 


$1 


1.005 (12) 


R,L,A 


1.008 (18) 


R,L,A 




1.26 (3) 


P 


1.152 (27) 


P 


$3 


1.63 (3) 




1.60 (5) 


C,B^ 


$4 


2.19 (8) 


A 


1.91 (5) 


Pd 


$5 


2.33 (8) 


C,B^,{T) 


2.19 (7) 


A 


$6 


2.25 (8) 


Pd 


2.29 (9) 


C,T 


2t 


M/gl 


Ops. 


M/gl 


Ops. 


$1 


1.266 (21) 


P 


1.143 (27) 


P 


$2 


1.274 (24) 


P 


1.161 (27) 


P 


$3 


2.32 (4) 


(P) 


2.38 (5) 


(P) 


$4 


2.43 (5) 


R,L,A 


2.51 (7) 


C,T 


^5 


2.61 (5) 


C,{T) 


2.45 (8) 


L,A 


$6 


2.78 (6) 


C,T 


2.72 (6) 


C 


$7 


3.23 (11) 


L 


3.17 (10) 


A 



Table 9: Mass estimates and dominant operator contributions in the Ol{l and 2^ chan- 
nels atT — 5Tc. The dominant operator types contributing are denoted with (without) 
parentheses if < (>)0.5. 
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T = 5T, 




= 9 


(3 = 12 


(3 = 1^ 


state 




L^-T = 243 


L^-T = 343 


L^-T = 483 




— 


— 


3.39 (15) 


— 


Ol 


1.352 (27) 


1.368 (28) 


1.342 (33) 


1.340 (30) 


Q-* 


2.34 (8) 


2.28 (6) 


2.33 (8) 


2.28 (5) 


n 


3.01 (37) 


3.65 (37) 


3.52 (37) 




K 


3.02 (26) 


3.16 (20) 


3.20 (20) 




It 


2.13 (3) 


2.04 (5) 


2.06 (6) 




ii 


2.02 (6) 


2.09 (5) 


2.06 (5) 




2+ 


2.69 (4) 


2.60 (10) 


2.64 (9) 




K 


2.52 (10) 


2.60 (9) 


2.54 (10) 





Table 10: Masses M/g^ for the other channels at T = bT^. 



an effect was encountered in the 1;^ sectors, as is apparent from comparing the two 
different parity states, which should be degenerate. This problem can only be solved 
with a larger basis including other operator types. For 2^ the basis equally consists of 
four operators of type A. In this case, however, we observe parity doubling with the 
ground state in the 2| channel (see below), and may thus be confident to have good 
projection on the 2^ ground state. 

At our largest lattice spacing, (3 = 9, we have performed an explicit check for finite 
volume effects. Comparison of the corresponding columns in Table |^ shows that the 
masses for those states that do not couple to winding operators (P, P^, T) are statis- 
tically compatible on ■ T = 243 L"^ - T = 3A^. It appears that already on the 
smaller volume most of our results for the spectrum are free of finite size effects, as 
was the case in the pure gauge theory. Of the operator eigenvectors projecting onto 
winding operators, those purely of type P or P^ may be discounted as mere finite vol- 
ume artefacts, they disappear from the spectrum in the infinite volume limit. More 
difficult are those states that couple both to non- winding, C or B±, and to winding 
operators, T, such as the 2^ states $5 and $6 whose masses are likely to suffer most 
from finite volume effects. At (3 = 9, where we have two volumes, we can resort to 
the larger lattice where the states are decoupled from the torelons, and trust that this 
number is as close to the infinite volume limit as the other states. On the finer lattices, 
however, we have only one volume, and hence it is likely that the corresponding masses 
are affected by the admixture of torelons to these states. 

As we have discussed, one important signal of finite volume and discretisation effects 
is the loss of parity doubling in the spectrum. On all our lattices we find within 
statistical errors a degeneracy of the and states not only for J = 1 (where it 
is expected even on the lattice) but also for J = 2. In the latter case, the relevant 
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0.05 0.1 

1/fi 



Figure 2: The approach to the continuum of the four lowest Ol[; states at T = ST^. 
Filled symbols denote glueballs, open symbols bound states of scalars. 

comparison is between non-torelonic states that do not couple to winding operators 
(and thus remain of finite mass in the infinite volume extrapolation). 

We have already commented that in extrapolating to the continuum limit there is a 
choice in the method of defining the lattice spacing. There is also some freedom in the 
functional form of the extrapolation. From the weak coupling expansion of the lattice 
action, the lowest order corrections to the continuum limit mass should be linear in 
the lattice spacing. Thus if we are sufficiently close to the continuum limit, a straight 
line fit to the data can be attempted. 

We have performed linear fits for the low lying and 2]^ states, using both tree 
level and non-perturbative definitions of the lattice spacing. All fits give acceptable 
X^/d.o.f . < 1. In Fig. ^ we show such fits using units of the bare coupling, and it is such 
a fit that we use for the ground state in the 2^ channel. Given the very large errors on 
the gradient for linear fits using the string tension to set the scale, we have insufficient 
data to decide if the latter, non-perturbative definition of the lattice spacing improves 
the scaling behaviour. 

The results for 0|; and 2^ are given in Table 0. The 2^ pure glueball couples to 
torelons on our finest lattice, so that its mass value there could not be used for an 
extrapolation. We quoted the value for (3 = 12 instead. We also compare with the 
corresponding glueball states in the pure gauge theory and find their masses to deviate 
at the percent level at most from those measured in the Higgs model. We conclude that 
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Olj; channel 


2'\_ channel 


gauge-Higgs 


pure gauge 


gauge-Higgs 


pure gauge 


scalar glueball 


glueball 


scalar glueball 


glueball 


1.002 (35) 

1.57 (9) 

2.22 (18) 

2.22 (25) 


1.58 (1) 
2.29 (3) 


2.46 (21) 

2.61 (5) 


2.62 (5) 



Table 11: Final mass estimates using the continuum extrapolation in the hare coupling 
of our data atT = bT^. Our data for the glueball are compared with the results obtained 
in the pure gauge theory 

we have a situation analogous to the fundamental Higgs model, with the pure gauge 
glueball spectrum being repeated, with additional bound states of adjoint scalars and 
little mixing between them. 

As comparison of Tables and |^ elucidates, the difference between the masses 
extrapolated to the continuum and those computed at finite lattice spacing are very 
small, of the order of a few percent for the lowest states, and mostly covered by the 
statistical errors of the values on the finest lattice. We therefore quote the latter for 
the other channels as well as our other parameter sets, noting, however, that for the 
heavier states this may well amount to an estimated 5% systematic error. 

4 Discussion 

With the numerical results for various correlation functions at hand, let us now discuss 
their implications for the spectrum of static correlation functions in four-dimensional 
SU(2) Yang-Mills theory at finite temperature. In order to express the masses entirely 
in units of temperature, we have to make use of the perturbative expression for the 
gauge coupling at that temperature scale, Eq. (^. Using this expression we have 
gi{2Tc) ^ 3.83r, ^^(2T,) ^ 3.07T and ^^(5TJ ^ 2.89T. The spectrum in these units 
is shown in Fig. ^. 

4.1 Comparison with four dimensions 

The first important question concerns the comparison with a four-dimensional simula- 
tion, in order to assess the quality of dimensional reduction in the temperature range 
under study. A detailed analysis of the low lying modes for SU(2) gauge theory has 
recently been presented for the temperature range 2Tc < T < ATc. Only the ground 
states in each quantum number channel are so far available in four dimensions, however. 



19 



15- 



10- 

M 
T 

5H 



0- 



T = 2Tr 



B 



0+070" 



(D CD 



T = 4Tr 



CD 4 



B 



B 



T = 5Tr 



a 



4 



4 



0) CD 



1+171+1" 



2+272+ 



0+070" 



1+171+1" 



2+272+ 



0+070" 



1+171+1" 



2+272+ 



Figure 3: T/ie spectrum of the lowest inverse static correlation lengths at each tem- 
perature. Filled symbols denote glueballs, open symbols bound states of scalars. 



To facilitate the comparison we briefly summarise the situation in four dimensions. 
There the lattice symmetry group is 1)4/1 ® Z{'^)rz where D^h is the point group of the 
discretised 2 + 1 dimensional space-time slice used to create screening state operators, 
and Rz is the Euclidian time reversal operator in the notation of [24|. It acts on the 
time component of the gauge field as Rz^q ~^o- Dimensional reduction converts 
the Aq component of the gauge field into an adjoint scalar field, and the scalar reflection 
symmetry is the "daughter" of Euclidian time reflection in the original theory. 

As the temperature is increased in the four-dimensional theory, the timelike direc- 
tion of the slice becomes squeezed, to the point where the slice becomes quasi-two- 
dimensional, and the effective symmetry group becomes C^. At this point the spectrum 
of screening masses should conform to a pattern predicted by this latter symmetry 
group, as has indeed been observed 0. If the perturbative matching of finite tempera- 
ture SU(2) pure gauge theory to the SU(2) adjoint Higgs model in three dimensions is 
correct, the screening masses should furthermore match those observed in this study. 
Hence we expect to see the 3 + 1 dimensional D^h states and A2 converging on the 
0"*" mass we measure (the quantum numbers Rz, R are suppressed here), and the 
and A2 on the 0~. The E states should agree with the J = 1. The and B2 should 
agree with the 2'*', and and B2 should have a common mass with the 2~. For a 
more detailed account of this symmetry breaking pattern we refer to 0. 

In P, |lOl and the above paragraph discussing the mapping of states, the behaviour 
under Euclidian time reflection Rz has not been specified. This leaves room for potential 
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ambiguities in comparing the states. For our comparison, we assume that the lowest 
states have been extracted in 0, without exphcit specification of their i^^-symmetry. 
Correspondingly, for the same spin and parity we compare with the lighter of our in- 
states. Further, because of parity doubling, has to be degenerate with Bf^ in the 



continuum limit. In [|10|, B^ still exhibits finite size effects, so that we take B^ instead 
to compare with our 2^ state. The four-dimensional data are taken from Table 4 in |^ 
and Table 2 in 0. 

The comparison of screening masses after these identifications is shown in Fig. 
Comparing the temperature dependence of the screening masses, we find agreement 
for the qualitative behaviour of the lowest states in all channels. Their temperature 
dependence is weak, with the 0~l slightly rising whereas all other ground states drop 
with T. In our simulation we find in addition that the temperature dependence becomes 
stronger for the higher excitations, cf. Fig. |^. But even on a quantitative level, we find 
a remarkable agreement for the low lying states. This could not have been anticipated 
given that the temperature is rather close to the phase transition, where we know the 
effective theory description has to break down. The only notable disagreement between 
the measurement in the full and the reduced theory is for the lightest spin 2 screening 
mass. But in comparing the heavier screening masses we have to remind ourselves 
that the effective theory is obtained by integrating out the non-zero Matsubara modes, 
whose smallest is of the order of ~ 27rT. The lowest spin 2 screening mass, on the 
other hand, is already larger than ~ 6T for the temperatures considered, and hence 
belongs to an energy region where the effective theory loses its validity. In this regime 
the non-local terms and terms of higher dimension which have been discarded in the 
construction of the effective theory are no longer suppressed, and disagreement for 
screening masses of that order and higher is in fact expected. 

We conclude that the three-dimensional SU(2) adjoint Higgs model in its metastable 
symmetric phase is the correct effective theory to describe the static correlation func- 
tions in a four-dimensional SU(2) pure gauge theory at finite temperature. The effective 
theory gives the correct qualitative picture, and for the screening masses in the expected 
range of validity of the effective theory even quantitative results at temperatures as 
low as T = 2Tc. 

4.2 The Debye mass 

Of particular interest for non-Abelian plasma physics is the Debye screening mass, 
which may be expanded as 



mo = + In ^ + CNgi + Oig-'T). 
The leading order result is known perturbatively |2^ to be 



^d"" = ^3 = ( y + ) gT. (19) 
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Figure 4: Comparison of screening masses as calculated in four dimensions directly 
(open symbols) and from the reduced three-dimensional theory used here (full symbols). 
In the point group notation of the open squares denote Af, diamonds Ai , circles 
Bi and triangles . 
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(l + 2)/<73' 


C2 


0{9'T)/gl 


T = 2T^ 
T = 4Tc 
T = 5T, 


1.12 (2) 
1.32 (2) 
1.34 (3) 


0.296 
0.360 
0.379 


1.14(4) 
1.14(4) 
1.14(4) 


-0.32(6) 
-0.18(6) 
-0.18(7) 



Table 12: The different contributions to the Debye mass, Eq. ^TB^), where (1+2) is the 
sum of the first two terms. 



At next-to leading order ~ g'^T, one can extract perturbatively a logarithm |]26|, but 
cn is entirely non-perturbative and has to be evaluated by lattice simulations. A gauge 
invariant, non-perturbative definition for m^) has been given in [El] and is based on 



Euclidian time refiection, R^. According to the definition in [24|, the Debye mass mo 



corresponds to the mass of the lightest state odd under this transformation, i.e. the 
lightest R = —1 state. According to our measurements this is the Ol ground state, as 
obtained from the operator B\\. As discussed in the previous section, our measurement 
is entirely consistent with the four- dimensional one, see Fig. ^. 

On the other hand, C2 has been determined directly in recent simulations [^, |2^ 
after integrating out Aq following a prescription in There the value C2 = 1.14(4) 
was found for the leading correction to the perturbative result^. 

Taking these data together, we are now in a position to assess the size of the various 
contributions in the expansion Eq. ([T8|) . These are summarised in Table 12. 



The emerging picture is quite interesting: the Debye mass is entirely non-perturb- 
ative, with its leading order value being much smaller than the g^T correction up to 
temperatures of the order of T ~ 10 ''Tc On the other hand, the 0{g^T) corrections 
are only 20% at temperatures as low as T = 2Tc, decreasing rapidly to 10% at T = 5Tc. 

We thus have to conclude that the scale dominating the physics of electric screening is 
the scale of the three-dimensional effective theory g^, or equivalently, the soft modes ~ 
g^T of the original finite temperature theory, contrary to the naive expectation ~ gT. 
The contributions on this scale are entirely non-perturbative. Contributions of higher 
powers in (yf, on the other hand, can be viewed as small perturbations. 



^ We remark that the operator i?| | has been simulated previously to determine the Debye mass . 
By fitting Eq. (^8|) to the temperature dependence of the mass of B\\, these authors find the term 
0{g'^T) to be consistent with zero, and upon subtracting the perturbative piece mj^^ determine the 
next-to-leading order correction to be C2 = 1.58(20). We consider this likely to be an overestimate, as 
it is two standard deviations above the results of . Adding the perturbative contribution and 

no other corrections, the resulting m^j likewise deviates from the four-dimensional result of and 
our present measurements. 
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5 Conclusions 



We have performed a detailed simulation of the spectrum of the SU(2) adjoint Higgs 
model in 2 + 1 dimensions at various points in its metastable confinement phase. The 
choice of the parameter values was motivated by the connection with four-dimensional 
SU(2) pure gauge theory at finite temperature, for which the adjoint Higgs model 
emerges as an effective theory for the static modes after perturbative dimensional re- 
duction. Our results are relevant both for the study of confinement in three-dimensional 
gauge theories, as well as for the question whether dimensional reduction is applicable 
to SU(2) pure gauge theory, and eventually to QCD. 

Regarding three-dimensional gauge theories, we found a dense spectrum of bound 
states, which consists of an almost unchanged replica of the glueball spectrum of pure 
gauge theory, to which additional bound states of adjoint scalars are added. The glue- 
ball and scalar states show very little mixing, and correspondingly the string tension 
is very close to its pure gauge value. The pure gauge quantities are furthermore in- 
sensitive to variations of the scalar parameters in the action. All these findings are 



completely analogous to the situation in the SU(2) fundamental Higgs model |]T5|. We 



conclude that the approximate decoupling of the pure gauge sector from the scalar 
sector is a very robust phenomenon in three dimensions, for which an explanation is 
lacking at present. 

Regarding its role as effective theory, we find that dimensional reduction works re- 
markably well down to temperatures of 2Tc, where the lowest screening masses com- 
puted in the effective and the full theories agree quantitatively. This agreement is 
rather striking for a temperature so close to the transition. These findings confirm the 
SU(2) adjoint Higgs model as the correct effective high temperature theory to describe 
the thermodynamics of four- dimensional SU(2) pure gauge theory. In particular, our 
investigation has established unambiguously that the Aq degrees of freedom for tem- 
peratures not too far above the phase transition constitute the lightest state in the 
three-dimensional effective theory, and hence may not be integrated out, as already 
discussed in It furthermore implies that the naive hierarchy of scales gT <^ g^T 

does not translate into a corresponding separation of the dynamics of the Aq and Ai 
in this temperature range. 

We expect no qualitative changes when moving from SU(2) to SU(3). An analogous 
investigation for SU(3) and three flavours of quarks poses no additional complications 
and appears desirable in the light of our results. 
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